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SUMMARY 

A  derivation  is  presented  of  a  point-source  method  based  on 
the  linearized  theory  for  obtaining  the  pressure  coefficient  on 
thin  wings  at  supersonic  speeds.  The  method  ia  applied  to  calcu¬ 
late  the  lift  distribution  of  a  thin  flat -plate  wing  having  a 
straight  leading  edge  and  an  arbitrarily  curved  wing  tip  whose 
elements  are  swept  behind  the  Mach  angle.  A  qualitative  basis  for 
obtaining  solutions  that  satisfy  the  Kutta-Joukowski  condition  is 
included . 

The  analysis  is  continued  to  formulate  the  velocity  potential 
for  regions  influenced  by  so-called  subsonic  trailing  edges  (that 
is,  edges  where  the  component  of  the  free-stream  velocity  normal 
to  the  trailing  edge  is  subsonic).  The  derivations  include  the 
solution  that  satisfies  the  Kutta-Joukowski  condition  or  any  of 
the  infinite  number  of  transition  solutions  involving  a  disconti¬ 
nuity  in  the  cross  velocity  (sidewash)  in  the  wake  of  a  subsonic 
trailing  edge.  As  an  example,  the  two  perturbation-velocity 
components  in  the  plane  of  a  trapezoidal  wing  and  the  upwash  over 
the  wing  tip  are  evaluated.  By  means  of  series  expansions,  the 
mathematical  nature  of  the  lift  distribution  in  regions  influenced 
by  arbitrary  subsonic  trailing  edges  is  indicated  when  the  Kutta- 
Joukowski  or  some  other  condition  is  applied  to  make  the  solution 
unique . 


IRTRODUCTIOK 

The  problem  of  determining  the  lift  distribution  of  arbitrary 
thin  wings  at  supersonic  speeds  has  been  partly  solved  by  use  of 
the  linearized  theory  in  references  1  and  2.  In  reference  1,  the 
point-source  method  of  reference  3  is  extended  to  include  the 
effective  sources  generated  by  the  slopes  of  the  streamlines  in 
the  region  between  the  wing  boundaries  and  the  foremost  Mach  waves. 
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The  velocity  potential  and  hence  the  aerodynamic  coefficients  in  the  , 

vicinity  of  wing  tips  at  supersonic  speeds  can  then  he  calculated. 

An  explicit  solution  for  the  slopes  of  the  streamlines  in  the 
region  included  between  the  wing  boundaries  and  the  foremost  Mach 
waves  is  presented  in  reference  2.  This  solution  allows  the  methods 
of  reference  1  to  be  extended  to  include  the  calculation  of  the 
velocity  potential  at  points  on  the  wing  influenced  by  interacting 
external  flow  fields. 

In  subsonic  aerodynamic  theory,  an  infinite  number  of  solutions 
will  satisfy  the  boundary  conditions  for  a  given  wing.  Some 
criterion  is  then  required  to  choose  from  the  manifold  of  solutions 
the  one  that  is  experienced  in  practice.  A  criterion  that  yields 
theoretical  results  in  approximate  agreement  with  experiment  is  the 
Kutta-Joukowski  condition.  This  condition  requires  that  the  stream¬ 
lines  leave  the  trailing  edge  of  the  airfoil  smoothly  (that  is,  that 
the  three  components  of  the  perturbation  velocity  are  continuous 
across  the  trailing-edge  boundary) .  « 

When  the  flow  is  supersonic,  the  Kutta-Joukowski  condition 
need  not  be  fulfilled  at  the  trailing  edge  of  the  airfoil.  Local-  * 

ized  compression  and  expansion  waves  equalize  the  pressures  from 
the  bottom  and  top  surfaces  behind  the  trailing  edge.  Nevertheless, 
because  of  the  close  similarity  of  supersonic  and  subsonic  flow 
when  the  component  of  the  free-stream  velocity  normal  to  the 
trailing  edge  is  subsonic  (the  so-called  subsonic  trailing  edge), 
a  strong  opinion  among  aerodynamicists  is  that  the  Kutta-Joukowski 
condition  will  also  approximately  hold  at  supersonic  speeds. 

As  indicated  in  reference  2,  the  formulations  of  reference  1 
may  be  extended  to  include  solutions  that  satisfy  the  Kutta- 
Joukowski  condition.  The  extension  shows  that  in  supersonic-flow 
theory,  as  in  subsonic-flow  theory,  an  infinity  of  solutions  can 
exist  for  a  given  wing  boundary.  Some  additional  criterion  is 
then  required.  For  subsonic  leading  edges,  the  added  boundary 
condition  applied  in  reference  1  is  that  the  perturbation-velocity 
components  shall  be  continuous  in  the  region  between  the  wing 
boundary  and  the  foremost  Mach  waves.  The  solutions  so  obtained  do 
not  satisfy  the  Kutta-Joukowski  condition  for  regions  influenced 
by  subsonic  trailing  edges.  If  a  discontinuity  in  the  cross¬ 
velocity  component  (sidewash  in  the  plane  of  the  wing)  is  per¬ 
mitted  (corresponding  to  the  shedding  of  a  vortex  sheet),  solu¬ 
tions  may  be  obtained  that  satisfy  the  Kutta-Joukowski  condition 
or  any  transition  solution  depending  upon  the  assumed  strength  of 
the  discontinuity  in  the  cross  velocity.  The  particular  g 

criterion  to  be  used  can  be  determined  only  by  experiment. 
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A  derivation  "based  on  linearized  theory  is  presented  in  the 
first  part  of  this  report  to  show  that  the  theoretical  pressure  dis¬ 
tribution  on  thin  wings  at  supersonic  speeds  may  be  obtained  from  a 
line  integral  and  a  surface  integral.  (A  line  source  was  applied  in 
reference  4  to  obtain  the  pressure  distribution  at  zero  lift.)  For 
thin  flat-plate  wings  or  for  the  lift  distribution  on  finite- 
thickness  wings,  only  line  integrals  are  required.  As  an  example, 
the  lift  distribution  is  evaluated  for  a  wing  having  a  straight 
leading  edge  intersected  by  an  arbitrarily  shaped  curve  whose 
boundary- line  elements  are  swept  behind  the  Mach  line  from  the 
wing-tip  and  leading-edge  intersection.  The  second  part  of  the 
report  formulates  the  expressions  for  the  velocity  potential  in 
regions  influenced  by  subsonic  trailing  edges. 

Examples  of  these  methods  are  presented.  The  solution  satis¬ 
fying  the  Kutta-Joukowski  condition  is  derived  for  a  swept  trap¬ 
ezoidal  wing.  Sidewash  and  upwash  velocities  in  the  vicinity  of 
the  wing  tip  are  also  calculated.  As  a  further  application,  by 
means  of  series  expansions,  the  mathematical  form  of  the  lift  dis¬ 
tribution  associated  with  the  wake  behind  a  general  subsonic 
trailing  edge  is  indicated.  When  this  result  is  combined  with  the 
line  integrals  in  the  first  part  of  the  report,  expressions  are 
obtained  for  the  lift  distribution  of  arbitrary-plan-boundary 
single  wing  tips  at  supersonic  speeds.  The  Kutta-Joukowski  and 
other  conditions  are  applied  to  formulate  feasible  solutions  for 
the  various  wing  areas.  The  research  was  conducted  at  the  NACA 
Cleveland  laboratory  during  August  and  September  1947. 


ANALYSIS 

Point-Source  Method  for  Evaluating 
Pressure  Coefficient 

A  convenient  approach  toward  evaluating  the  pressure  coeffi¬ 
cient  on  thin  finite  wings  at  supersonic  speeds  is  to  establish 
the  velocity  potential.  A  point-source  integral  expression  for 
the  velocity  potential  at  any  point  on  thin  wings  at  supersonic 
speeds  is  given  in  references  1  and  3  as 
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where  here 


cp  perturbation-velocity  potential  on  top  wing  surface 

U  free-stream  velocity  considered  parallel  to  x  or  i  axis 

O  slope  of  streamlines  near  x,  y  plane  measured  in  y  =  constant 

planes 

t,  I  Cartesian  variables  of  integration,  x  and  y  directions, 
respectively 

x,y  Cartesian  coordinates  of  point  where  cp  is  evaluated 

3  =  ^|m2-i 

M  free-stream  Mach  number 

(A  complete  list  of  symbols  used  in  this  report  is  presented  in 
appendix  A.)  In  equation  (l),  the  integration  of  a  is  conducted 
in  the  x,y  plane  for  the  disturbed  flow  field  S  included  in  thG 
forward  Mach  cone  from  the  point  (x,y).  For  example,  the  integra¬ 
tion  for  the  wing  shown  in  figure  1(a)  is  over  the  area  Stf.  In 
regions  influenced  by  a  wing  tip,  the  disturbed  flow  area  between 
the  wing  boundary  and  the  foremost  Mach  line  must  also  be  included. 


If  the  velocity  potential  is  computed  at  a  point  (x+dx,y),  as 
in  figure  1(b),  the  area  enclosed  by  the  curve  efgc'  is  the  same 
as  the  area  abdc  and  the  relation  of  the  distances  in  the  denomi¬ 
nator  of  the  integrand  is  the  same  as  was  considered  at  point  (x,y). 

The  slopes,  however,  will  have  the  value  a  +  ^  dx.  In  addition, 

the  area  a'bd'gfe  will  contribute  to  the  velocity  potential.  The 
value  of  cp  then  becomes 


cp  +  ^3?  dx 


Cr 


JJ 


(a+||  dx)d|dii 
/\j(x-i)2-P2(y-r])2 


El+dx 


adndl 


JJ 


|x  ^U-!)2-32(y-T)) 


(2) 


where  is  the  value  of  £  along  the  leading  edge.  Integration 

yields 
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cp  + 


Sep 

s 


dx 


nr 


^ J 


^(x-U2-P2(y-n)2 


adTjdx 


Jabd  rf(x-i)2-P2(y-r})2 


(2a) 


where  abd  represents  the  path  of  integration.  Subtraction  of 
equation  (l)  from  equation  (2a)  gives 


* 


r 


%>  =  m  U 
)x  Jt 


rr 


j§^ 


gdri 


abd  i^(x-|  )2-P2(y-T)) 


(3) 


where  |  and  (7  are  evaluated  in  the  line  integral  as  functions 
of  t]  along  the  line  ahd  of  figure  1(a). 


The  x  component  of  the  perturbation  velocity  may  then  be 
computed  as  the  surface  integral  of  the  rate  of  change  of  the  stream¬ 
line  slopes  with  respect  to  i  plus  a  line  integral  along  the 
foremost  boundary  of  the  disturbed  flow  field  included  in  the 
forward  Mach  cone  from  the  point  (x,y) .  The  surface  integral,  of 
course,  vanishes  in  regions  where  O  is  constant.  (Evaluation  of 

the  surface  integral  along  lines  of  discontinuous  ^  results  in 
additional  line  integrals.) 

For  some  problems,  the  integrations  of  equation  (3)  may  be 
more  easily  evaluated  in  a  set  of  oblique  coordinates  (u,v)  whose 
axes  lie  parallel  to  the  Mach  lines.  The  transformation  equations 
are 
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u  = 

v  =  §(U  Pn) 

1 

II 

T1  =  i(v-u) 

°w  =  ^(i-  Py) 

vw  =  ^(2+Py) 

x  =  m(tv+«w) 

y  =  s(tw_uw) 

(4) 


Inasmuch  as  the  elemental  area  in  the  (u,v)  coordinate  system 
is  ^fdudv,  equation  (3)  may  be  written  as 


M 


&  =  _  JL  ff  +  ^)dudv 
2P«  JJ  g  Y (uw-u) (vw-v) 


U 

2  Pit 


g  (dv-du) 


abd  ^(uw-u)(vw-v) 


(3a) 


where  (uw,vw)  corresponds  to  the  point  (x,y)  and  the  evaluation 
of  the  line  integral  is  again  conducted  along  the  foremost  line  of 
disturbances  included  in  the  forward  Mach  cone. 


When  only  one  wing  tip  Is  included  in  the  forward  Mach  cone 
from  the  point,  the  surface  integral  over  the  area  Sp  (fig.  2(a)) 
generated  by  the  upwash  over  the  wing  tip  may  be  replaced  by  an 
equivalent  integration  over  a  portion  of  the  wing.  This  simplifi¬ 
cation,  which  follows  from  the  derivations  of  reference  1, 
eliminates  the  necessity  of  evaluating  the  slopes  of  the  streamlines 
in  the  external  flow  field  Sp.  The  velocity  potential  at  point 
(uw,vw)  on  the  wing  of  figure  2(a)  is  obtained  from  reference  1  as 


0  ydudv 

VK-)  (TW~V) 

(ag-KTipJdudv 
2  /^(uw-u)(vw-v) 


* 


* 


# 


(5) 
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where  and  Qiji  are  the  bottom-  and  top-surface  wing  slopes 9 

respectively. 


Because  evaluation  of  cp  from  equation  (5)  requires  the  applica¬ 
tion  of  separate  sets  of  wing  slopes  in  the  respective  regions 
and  the  calculation  of  ^  will  result  in  a  line  integral 

along  the  dividing  boundary.  The  derivation  presented  in  appendix  B 
then  yields  equation  (6). 


m 


uw-u)(vw-v) 


rvN 


u 

4  Pit 


®w,2 


5(c?g+ap)  d(a-g-K7ip) 
5u  +  dv 

^(uw-u)(vw-v) 


dudv 


U 

2  (Bn 


U 

2p7T 


ab 


bd 


aT(dv-du)  U 

Y(uw-u)(vw-v)  2Pn 


bod 


l- 


du2(yw) 


dv 


w  j 


dv 


2  ][( uw-u)(vw-v) 


(aB-K7T)  (dv-du) 

2  ty(uw-u)(vw-v) 


(6) 


du2(vw) 

wnere  — - - 

dvw 

at  v  =  vw. 


is  the  derivative  of  the  wing-houndary  curve  u  =  U2(v) 


In  order  to  obtain  the  lift  distribution  as  a  function  of  angle 
of  attack  a,  only  the  solution  for  a  thin  flat  plate  is  required. 
Because  C£  =  -GT  =  a,  only  the  two  line  integrals  along  the  curves 
ab  and  bd  remain.  The  portion  of  the  equation  (6)  that  yields  the 
lift  distribution  is  then 
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2P« 


U 

2P« 


-a  (dv-du) 


Jab  V(aw'u)(vw“v) 


bd 


(vw-v) 


(7) 


If  the  equations  of  the  wing  boundaries  (fig.  2(b))  on  the  two 
branches  separated  by  the  origin  are  denoted  by 

y  =  v1(u)  or  u  =  u^(v) 

and 

v  =  V2(u)  or  u  =  u2(v) 

equation  (7)  may  be  partly  evaluated  to  yield 


n 


dcp  _  Ua 
^x  “  2fbt 


(dv-du) 

ab 


.  Ua 

„  d^2(vw) 

A  vw"vl(u2(vw)) 

1  dvw 

h 

DO 

? 

(7a) 


where  v^  is  a  function  of  u2(vw).  The  expression  for  ^ 

represented  by  equation  (7a)  is  applied  in  the  second  part  of  this 
report  to  derive  solutions  for  the  lift  distribution  that  satisfy 
the  Kutta-Joukowski  condition  for  general-plan-fonn  single-wing- 
tip  boundaries. 

As  an  illustration  of  equation  (7a),  the  pressure  coefficient 
on  the  top  wing  surface  as  given  by  the  relation 


CP  =  -^ 


(8) 


will  be  calculated  for  the  wing  shown  in  figure  3 
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Inasmuch  as  the  leading  edge  and  the  wing  tip  satisfy  the  equa¬ 
tions  v  =  -  kqu  and  u  =  Ug(v),  the  integral  of  equation  (7a) 
becomes 


For  comparison,  equation  (ll)  is  also  derived  in  appendix  C  as  equa¬ 
tion  (C4)  by  the  method  of  reference  1. 


Equation  (ll)  shows  that  the  pressure  coefficient  assumes  the 
Ackeret  value  for  a  swept  wing  along  the  innermost  Mach  line  from 
the  origin  (vw  =  0)  and  approaches  infinity  along  the  subsonic  leading 
edge  of  the  wing  tip.  The  portion  of  the  solution  that  gives  the 
infinity  (that  is,  equation  (10))  is  generated  by  the  line  integral 
along  the  Mach  line  bd  in  equation  (7).  This  quantity  is  zero  when 

dupCv-^) 

the  sweep  on  the  wing  tip  is  77  or  — — -  =  1-  A  change  in  sign 

u  nvTj 

du.2  w  N 

occurs  when  - — >  1  (regions  influenced  by  subsonic  trailing  edges) 
dvw 

in  equation  (10),  which  leads  to  a  reversal  of  sign  in  the  pressure 
coefficient  over  a  portion  of  the  top  wing  surface.  This  pressure 
reversal  leads  to  lost  lift  and,  because  of  trailing-edge  suction,  to 
decreased  lift-drag  ratio  as  compared  with  the  solution  that  satisfies 
the  Kutta-Joukowski  condition. 


The  pressure  coefficient  that  satisfies  the  Kutta-Joukowski  con¬ 
dition  along  the  subsonic  trailing  edge  must  match  the  solution  given 
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In  equation  (ll)  along  the  constant-vw  line  evaluated  by  the 
du2 

relation  =  1.  The  pressure  coefficient  must  furthermore  be 

zero  along  the  subsonic  trailing  edge.  A  function  that  meets  both 
of  these  requirements  is 


2a  (kl+1) 

Hfil 


tan“^- 


kl[uM-u2(vy)] 

vw+kiu2(vw) 


(12) 


The  wake  from  the  subsonic  trailing  edge  alters  the  lift  distribu¬ 
tion  over  a  portion  of  the  wing  in  such  a  manner  that  the  discon¬ 
tinuous  term  in  equation  (ll)  is  canceled.  In  the  second  part  of 
this  report,  equation  (12)  is  shown  to  be  the  correct  solution  for 
the  pressure  coefficient  in  regions  influenced  by  subsonic  trailing 
edges  of  a  wing  whose  leading  edge  is  straight  (fig.  3)  when  the 
Kutta-Joukowski  condition  is  to  apply. 


Formulation  of  Velocity  Potential  Including  Effects 
Associated  with  Wake  behind 
Subsonic  Trailing  Edge 

In  order  to  clarify  the  present  discussion,  a  brief  review  of 
the  concepts  presented  in  reference  1  is  included.  The  calculation 
of  the  velocity  potential  at  some  point  on  the  surface  of  thin  wings 
at  supersonic  speeds  requires  a  knowledge  of  the  streamline  slopes 
(a  of  equation  (l))  included  in  the  forward  Mach  cone  from  the 
point.  The  slopes  of  the  streamlines  over  the  area  Sw(i+2) 

(fig.  2(a))  are  just  the  wing-section  slopes.  The  slopes  of  the 
streamlines  in  the  region  Sj)  must  then  be  calculated,  either 
implicitly  or  explicitly. 

The  slopes  of  the  streamlines  in  the  region  Sj)  may  be  evaluated 
by  assuming  that  a  thin  impermeable  diaphragm  coincides  with  a  stream 
sheet  in  the  plane  of  the  wing.  This  assumption  does  not  alter  the 
flow  over  the  wing.  The  diaphragm  may  then  be  regarded  as  an  exten¬ 
sion  of  the  wing  to  isolate  the  top  and  bottom  surfaces. 

The  coincidence  of  the  diaphragm  with  a  stream  sheet  requires 
that  no  pressure  difference  exist  between  its  top  and  bottom 
surfaces.  This  condition  may  be  satisfied  by  requiring  either  that 
the  surface  velocity  potential  be  continuous  in  passing  from  the  top 


886 


MCA  TN  No.  1585 


11 


to  the  bottom  surface  of  the  diaphragm,  as  in  reference  1,  or  that 
the  component  of  the  perturbation  velocity  along  the  free -stream 
direction  (the  x  component)  be  continuous. 

If  the  x  component  of  the  perturbation  velocity  (obtained  by 
differentiating  the  velocity  potential)  is  equated  on  the  top  and 
bottom  surfaces  of  the  diaphragm,  the  defining  equation  for  the 
slope  A  of  the  stream  sheet  SB  (fig.  4(a))  is  obtained  as 


* 


s 


rr 


Sd 


^(l,T])d|dq  a 

r* 

(aB-oT)d|dq 

ty(x-S)2-P2(y-n)2 

^  2  A/(x-e)2-82(y-q)2 

Integration  of  equation  (13)  yields 


(13) 


SD 


V(x-l)2-P2(y-n)2 


rr 


w 


(Ob-°T)  didr] 

^(x-£)2-P2(y-Tj)2 


(14) 


If  the  function  f(y)  that  results  from  the  partial  integration  is 
zero,  the  equation  derived  in  reference  1  for  defining  X  is 
obtained . 


A  calculation  of  the  velocity  potential  on  the  top  surface  of 
the  diaphragm  for  any  value  of  the  function  f(y)  may  be  made. 
Again  with  reference  to  figure  4(a),  the  velocity  potential  becomes 


CTd|dt} 


JJ 

rr 


Sw  ^(x-  l)2-P2(y-n) 

Xd|dti 


’D 


f\j(x-  D2-P2(y-TiT 


Substitution  of  equation  (14)  into  equation  (15)  yields 


(15) 


rr 


V 

*  =  "  * 


(oB-KJT  )d|dq 


+  f  (y) 


(16) 


Vu-D2-P2(y-n)2 
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A  similar  calculation  of  the  potential  on  the  bottom  surface  of  the 
diaphragm  changes  only  the  sign  affixed  to  f(y).  The  function  f(y) 
thus  represents  one -half  the  difference  in  potential  between  the  top 
and  bottom  surfaces  of  the  diaphragm.  In  the  case  of  thin  plate 
wings ,  Og-KTiji  =  0  and  f(y)  is  the  potential  on  the  top  surface  of 
the  diaphragm. 

Despite  a  discontinuity  in  the  velocity  potential  across  the 
diaphragm,  the  x  and  z  components  of  the  perturbation  velocity 
are  continuous.  This  apparently  contradictory  state  of  affairs  may 
be  rationalized  by  the  presence  of  a  vortex  sheet  in  the  plane  of 
the  wing.  The  strength  of  the  vortex  sheet  may  be  chosen  (by 
controlling  f(y))  to  yield  solutions  for  the  velocity  potential 
that  satisfy  the  Kutta-Joukowski  condition  for  subsonic  trailing 
edges. 


The  function  f(y)  is  related  to  the  strength  of  the  vortex 
sheet.  If  no  wing  nor  material  object  is  intercepted  (for  the 
values  of  y  under  consideration)  to  generate  a  wake,  the  strength 
of  f  (y)  remains  unaltered  for  all  values  of  x. 

In  the  equations  that  follow,  two  simplifications  are  introduced 
(l)  all  the  derivations  are  conducted  in  the  oblique  coordinate 
system,  and  (2)  the  derivations  consider  only  thin  flat-plate  wings. 
Neither  simplification  decreases  the  generality  of  the  solution. 

The  results  for  the  actual  wing  at  angle  of  attack  may  be  evaluated 
as  the  sum  of  the  results  at  zero  angle  of  attack  and  the  solution 
at  angle  of  attack  for  the  thin  flat  plate  of  the  same  plan  form. 

The  external  flow  field  may  be  divided  Into  two  parts  (as  in 
fig.  4(b)),  separated  by  the  value  of  y  that  denotes  the  Junction 
between  the  leading  and  trailing  edge.  The  foremost  Mach  line 
originating  on  the  leading  edge  generally  represents  a  line  of 
infinitesimal  disturbance  along  which  f (y)  may  be  set  equal  to 
zero.  The  function  f(y)  will  remain  zero  for  values  of  x  that 
do  not  intercept  the  wing.  The  velocity  potential  (in  the 
z  =  0  plane)  is  thus  zero  in  the  field  The  Sj)  i,  on 

the  other  hand,  lies  in  the  wake  of  the  wing.  In  this  region,  f(y) 
will  generally  be  other  than  zero  (corresponding  to  the  generation 
of  a  vortex  sheet)  and  may  be  evaluated  along  the  wing  trailing 
edge. 


If  the  region  Sp  i  Is  temporarily  considered  as  a  portion  of 
the  wing,  the  method  of  reference  1  may  be  applied  to  evaluate  the 
potential  on  the  top  surface  of  the  diaphragm  or  on  the  wing 

surface.  The  virtual  wing  tip  is  then  the  y  =  constant  line 
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denoting  the  junction  between  the  leading  and  trailing  edges.  Inas¬ 
much  as  f(y)  is  zero  in  the  region  "the  effects  of  this 

disturbed  flow  field  may  be  evaluated  by* the  methods  of  reference  lj 
equation  (5)  applied  to  the  virtual  wing  boundary  will  result.  When 
the  flow  fields  are  subdivided  and  renumbered,  as  on  figure  4(c), 
the  contribution  to  the  velocity  potential  of  the  disturbed 
fields  Sj)( 2+4+3)  and  S^2  ls  given  by  equation  (5)  as 


(oB-KJT)dudv 

ty(iVu)(vw-v) 


(X3-^3)dudv 

l/(uv-u)(vw-v) 


because  oB  =  -Cip.  Only  the  slopes  of  the  streamlines  in  the 
regions  S^q  and  Sq^q  need  then  be  considered.  The  velocity 
potential  thus  becomes 


The  origin  of  the  coordinates  for  the  wing  of  figure  4(d)  is 
placed  at  the  junction  of  the  subsonic  leading  and  trailing  edges. 
The  wing  boundaries  are  defined  on  the  two  sides  of  the  origin  by 
the  two  sets  of  equations: 

v  =  vq(u)  or  u  =  uq(v) 

v  =  vg(u)  or  u  =  u2(v) 

dug 

(It  will  also  be  implicitly  assumed  that  >  1.)  Insertion  of 

the  limits  of  integration  into  equation  (17)  yields 


(18) 


14 


MCA  TN  No.  1585 


By  use  of  Abel's  solution  (reference  2,  5,  or  6),  equation  (18)  may 
be  inverted  to  yield 


rv2(u) 

ryp 

-adv 

Xdv 

M  d 

T 

Vl  (u)  ^Jvp-v  w 

^2(u)  ^fvp-v 

U  3u 

•J 

f  (y)aup 

p; 


(19) 


The  velocity  potential  at  points  on  the  wing  influenced  by  the 
subsonic  trailing  edge  (fig.  4(e))  is 

rrv2(u) 


ru2(vw) 


_  U 

^  ~  ~  Mit 


Mjt 


du 


pv 


u2(vw) 


uw-u 


du 


-ccdv 

LJvi(u)  ^ 
r* w 


Xdv 


Jv2( 


„)  ifl 


vL(u) 


-adv 

») 


(20) 


The  second  member  of  equation  (19)  may  replace  the  first  member 
along  lines  of  constant  vp  (or  vw)  that  extend  across  the  wing. 
Equation  (20)  then  becomes 


ru2(vw) 


-'Vw 


(21) 


Equation  (21)  gives  the  velocity  potential  on  the  wing  in  regions 
influenced  by  subsonic  trailing  edges.  Inasmuch  as 

(Vp-UjA 

— — ]  may  be  arbitrarily  chosen,  it  is  apparent  that  an 

infinite  number  of  solutions  can  satisfy  the  boundary  conditions. 
If  f  =  0,  equation  (21)  reduces  to  the  result  that  would  be 
obtained  by  the  methods  of  reference  1.  On  the  other  hand,  if  the 
Kutta-Joukowski  condition  is  to  apply,  the  function  f  may  be  so 
chosen  that  the  flow  leaves  the  subsonic  trailing  edge  smoothly. 


After  the  function  f  has  been  chosen,  the  second  member  of 
equation  (19)  may  be  evaluated.  The  slopes  of  the  streamlines  X  may 
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then  be  calculated  by  use  of  Abel's  equation.  These  slopes  in 
general  consist  of  two  parts,  A  =  Aa+A^>  which  satisfy  the  two 
equations 


rv2(u) 


-adv 


Jv 


l(„)  f**  Jv2(u) 


D  A^dv 


=  0 


and 


r*  v 


^  Aadv 


V2( 


»)f^ 


M  c> 
U  3u 


f  (y)duj) 

,D  f  ^ 


Equation  (22)  has  been  evaluated  in  reference  2  to  give 


2a 

A 

v2(u)-vx(u)  _  _n  I\ 

v2(u)-vi(u) 

it 

L‘ 

v-v2(u)  tan  ^ 

v-v2(u) 

(22) 


(23) 


(22a) 


Equation  (23)  may  be  solved  by  Abel's  equation  (reference  2,  5,  or  6) 
to  give  r  , 

/vp-uD\ 

V  M  ]' 


=  - 


_M_  c> 
Uit  c5v 


'-u 

pv 

a 

L  J 

VD 

. 

v2(u) 

|dup 


dv- 


P 


(23a) 


In  this  manner,  the  slopes  of  the  streamlines  in  the  region  Sp^ 
of  figure  4(b)  may  be  evaluated.  If  the  region  Sp  ^  is  treated 
as  a  portion  of  the  wing,  the  slopes  of  the  streamlines  in  the 
region  Sp^2  be  evaluated  from  the  equations  of  reference  2. 


Examples  of  Solutions  That  Satisfy 
Kutta-Joukowski  Condition 

Swept  trapezoidal  wing.  -  As  a  simple  example,  the  theory  of 
the  trapezoidal  wing  shown  in  figure  5  will  be  outlined.  The  leading 
and  trailing  edges  of  the  tip  are  defined  by  the  equations 
v  =  V]_(u)  =  -  k]_u  and  v  =  vgu  =  kgu,  where  ki>  0  and  0  <  k2  <  1 
Because  the  wing  boundaries  are  straight  lines,  conical  flow  will  be 
assumed . 
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Because  conical  flow  has  Been  assumed  in  the  region  influenced 

by  the  wing  tip,  the  cross  velocity  -r-  must  be  either  constant  or 
/  °y  Sf  df 

a  function  of  y/x.  Because  f  is  a  function  only  of  y,  ^ 
is  constant  in  the  region  Sj)  l.  The  function  f(y)  is  then  of  the 
form  &q  +  a]_y  where  aQ  and  a^  are  constants.  Inasmuch  as  f(y) 
must  be  zero  along  the  line  y  =  0,  ao  is  zero.  Equation  (19)  then 
becomes 


■odv 


■*!<! 


r’D 


Jk2u 


A  dv 

ftp 


M  3 

u  ST 


al(vD-uI>)duD  2ax  - 

=  —  ^|U-VB  C24' 


VD 


M 


^|a-uj) 


Equation  (21)  then  yields  the  velocity  potential  as 


/~v. 


2a^ 

V  =  "  msc 


w/k2 


■  U. 


Jy ; 


w 


^u-vw  ,  U 
a'm 


V 


r*v. 


vw/k2 


du 

^uw-u  J 


v 


-k^u 


-ccdv 


(25) 


The  two  integrations  indicated  in  equation  (25)  may  be  evaluated  by 
formulas  111  and  113  of  reference  7  to  give 

rvw/k2 


CPfl.  = 


2s*j 

Mrt 


/Ju- 


'w 


V'v-* 


du  = 


2a^ 

Mit 


/tfl-k2  , - 

^vw(k2uw-vw) 


+  (vw-uw)  tan-1 


Vw(1-k2) 

k2uw-vw 


(26) 


and 


ny 


du 


,w/k2  ^J.klu  fP* 


w 


ccdv 


2Ua 

Mn 


A  kn  +ko  I  ~  ; 


klUW+VW  J. _ -1 

+  -  tan  x 


fki 


^(kgUy-Vj 
I  vw(kx+k2) 


(27) 
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The  potential  cp  =  cpa  +  cp^  must  assume  the  value  aqy  along  the 
curve  vw  =  k2uw.  Equation  (25)  (from  equations  (26)  and  (27)) 

meets  this  requirement.  Moreover,  must  he  equal  to  aq  aloi 

the  curve  vw  =  k2uw  if  the  Kutta-Joukowski  condition  is  to  he 
satisfied.  Differentiation  of  equations  (26)  and  (27)  yields 


W 


fa 

Jt 


l-k. 


1 


'v 


k2Vvv 


(k2+l)  +  2tan' 


(28) 


_  Ua 

(k2+l)  kq+k2  k 

.... v»  ' ,  a-6!1  « 

— 

kl(k2%-vw) 

Sy  it 

*2  \ 

vw(ki+k2) 

„  dcp  ^a  ^Pb 

Setting  =  -gj-  +  =  aq  at  vw  =  k2u„.  gives 


(29) 


. I  kl+k2 

ai--uaAj^4 

The  cross  velocity  on  the  wing  is  then 

(l-kq) 


Ua 

n 


tan 


-1 


1 


kl(k2uv-vv) 

Vw(kq+k2) 


(30) 


-  2 


kl+k2  x  _q 
— r-r —  tan 
l-k2 


In  a  similar  manner, 

S' 


dcp  aCpa  ^h  Ua  (kl+l)  _  _q  Af^ ^V1^) 
^=^T  +  lF  =  ^-^-tan  l]j  vw(kq+k2) 


(32) 


which  is  in  agreement  with  equation  (12)  and  with  the  results  of 
reference  8.  W.  D.  Hayes  of  North  American  Aviation,  Inc.,  has 
obtained  equation  (32)  from  conical-flow  relations. 

Both  the  x  and  y  components  of  the  perturbation  velocity 
are  continuous  (equations  (31)  and  (32))  across  the  subsonic  trailing 
edge,  as  required  by  the  Kutta-Joulcowski  condition.  The  upwash  over 
the  wing  tip  is  also  continuous.  From  equation  (22a), 
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■>  2a 

\  =  T 


D(k1+k2) 


.|up(k 

^  VD* 


'k2uP 

From  equations  (24)  and  (23a), 

X  _  2al  d 
Aa  W 


-  tan 


-1 


1 


uD(ki+k2) 


vP-k2uP 


(33) 


A 

|u-vp 

k2u 


^v-vD 


dvD 


Evaluation  of  Xa  and  replacement  of  u  and  v  by  up  and  vp, 
respectively,  gives 


^a(uD>TD)  =  "  if 


ki+k2 


1-k. 


<  TT 


up(l-k2) 


l 


k2uD 


+  log 


*JUD~VD 


jjvD-k2uD  +  jjuD(l-k2) 


+ 


Addition  of  A  a 


_ vp-up _ ' 

2  /jvp-kjjup  +  ^Jupi(l-k2)J  > 

and  Xq,  yields 


(34) 


2a 

+««-i  A 

(kq+k2)uD  * 

kl+k2  /^up-vp 

it 

TioJi  *  \ 

vp-k2up  1  \ 

2  ^Vp-kgUq,  +  uD(l-k2) 

(35) 


Equation  (35)  is  valid  only  in  the  region  Sp^q  of  figure  5.  As 
v  approaches  k2u,  X  from  equation  (35)  approaches  -a,  as 
required  by  the  Kutta-Joukowski  condition.  The  solution  then  satis¬ 
fies  the  Kutta-Joukowski  condition  in  all  the  perturbation-velocity 
components  even  though  only  one  of  the  components  was  applied  to 
evaluate  the  constant  aq. 


The  x  and  y  components  of  the  perturbation  velocity  and  the 
slopes  of  the  streamlines  in  the  region  Sp^q  are  plotted  as  a 
function  of  Py/x  in  figures  6  to  8  for  a  wing  on  which  kq  =  2 
and  k2  =  0.5.  The  x  and  y  components  (figs.  6  and  7)  of  the 
perturbation  velocity  decrease  as  Py/x  is  changed  from  -1  (inner¬ 
most  Mach  line)  to  -l/3  (the  wing  tip).  Both  the  x  and  y  com¬ 
ponents  remain  constant  in  the  region  “  \  -  if  “  0  •  The  streamline 
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slopes  (fig.  8)  ©re  equal  to  -a  on  the  wing.  The  slopes  change 

1  Ry 

from  -a  to  +»  in  the  region  -  —  ^  <  0.  Similar  values  were 

found  by  Lagerstrom  in  an  investigation  for  the  Douglas  Aircraft  Co. 

Generalized  wing  tip.  -  For  this  example,  a  function  f(y)  that 
can  be  expanded  in  a  power  series  may  be  chosen  for  the  region  Sj^  i 
of  figure  4(b).  The  function  f(y)  may  then  be  expressed  as 


Also,  from  formula  750  of  reference  7 , 


The  integrand  in  the  second  member  of  equation  (19)  then  becomes 
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Integration  of  equation  (38)  with  respect  to  uj)  yields 


P  £.(/.)dup  _  .  J  _ 

al  a0 

3 

(u-vd)2 

j  m 

JVj)  11 

- -  -  v_d; 

[u-uD 

M  2(u-vd) 

2 

5 

1 

a2  al  1X3  a0 

(u-vD)2 

T 

m2  2M(u-vd)  '  2X4  ^  .2 

3 

a3 

&2  1X3 

1x3x5  a0 

1 

CO 

M 

2M2(u-Vp)  2X4  M(u-v-q)2 

2x4x6 (u-vp)^ 

4 

(39) 


Equation  (39)  may  he  simplified  "by  collecting  the  coefficients  of 
the  constants  aQ,  a^,  ag,  etc.: 


ru 


VD 


f(y)duD  2  ( ,  1  .  1x3  1X3X5 

=  a0  (u-vD)6  1  +  + 


2X2  '  2X4X3  +  2X4X6X4  +  *  *  * 


1  / 

al  _  \2  /  1  .  1  1X3  .  1X3X5 

""  ^  (a"*yj))  (  O*  .*7  ^  *  +  f~i.  .  m  .  .n.  .i—  +  .  •  . 


2  2x3  2x4x4  2x4x6x5 


+  “§  (u-vD)2 


M 


M‘ 


A  _i_ 

^  3  +  2x4 


1X3  1X3X5 

2X4X5  +  2X4X6X6  + 


(u 


\z(l  1  1X3  1X3X5  \ 

D'  \4  2X5  2X4X6  +  2X4X6X7  +  *  *  *  ) 


+  •  •  • 


Now  by  formulas  750  and  482  of  reference  7 , 
n  1 

1  _ 1  1X3  1x3x5 

~  n  2(n+l)  +  2x4 (n+2)  +  2x4x6 (n+3)  +  * 

_n»)r]D 

r(n  +  5) 


0 


(40) 


(41) 
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Equation  (40)  then  becomes 
nx 

f(y)duD 


ao  (u-vd)2xT“1T(u-vd)2x31: 


YD 


2X2 


a2 


222  2X2X2X2 


a3 


■I 

m=n 


(-1) 


m 


fm 


2m+l 

(u-vD)  2 


1~(^1)r  (l) 

rHi) 


Partial  differentiation  of  equation  (42)  with  respect  to 


ru 


f(y)duD 


a0 


1 

(u-vj))2 


fa 

M 


i 

v  \2  1.  f2  ,  .2  2X1 

d)  IP  (u-vi>)  jt 


5 

/■„  „  \2  3X2X1 

(  Td)  m  + 
2  2X2 


oo 


m=0 


2m- 1 

(-l)m  ^  (u-vD)  2 

VP 


r^)  r(l) 

TF7W 


Equation  (43)  shows  that  either  member  of  equation  (19)  is 
only  of  the  quantity  u-vD. 


P 


(42) 

yields 


(43) 

function 
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The  portion  of  the  velocity  potential  associated  with  the 
function  f(y)  Is  given  from  equation  (21)  as 


where  v^-  has  replaced  vp  inside  the  bracket,  because  only  u  and 
uD  are  variables  in  the  manipulation.  From  formula  750  of  refer- 


Substituting  equations  (43)  and  (45)  into  (44),  with  vw  replacing 
Vp  and  U2(VW)  shortened  to  U2,  yields 
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=  2(u2-vw)2 


a0 


(uw-vw)2 


(u2-vw)2 


2a 


3  1 

2(uw-vw)2  M(uw-vw)2 


(u2-vw) 


1X3  _ ^0 

2X4 


1X2  &i 


a2 


2X1 


5  3  M  13  1 

(uw“vw)2  2(uv“Vv)2  M2(uw-yv)2  2s 


(u2-vw)2 


1X3X5 


1X3X2  al 


2X4X6  2  5  M 

(uw-vw)2  2x4(uw-vw)2 


1X2X1  a2  3X2X1  a3 


—  M2  — X— X-  M3  — 

31  2  2X2X2  (n-v)2 

2>^x|(uw-vw)  V 


(u2-vw)2 


1X3X5X7 


Ix3x5x2  al 


2x4x6x8  9  5  M 

(uw-vw)2  2x4x6 (uw-vw) 2 


(46) 


Collecting  the  coefficients  of  the  constants  slq,  a^,  a2 
gives 
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2x2x1  a2  (u2"vv)2 


3xi 

2  2 


M‘ 


(uw-vw)2 


2X3X2X1  a3  (u2~Tw)2 
,3 


5X3X1  M' 
2  2  2 


1 

(uw-vw)2 


.  1 
2X7 


(u2~Vy)  1X3  /u2~vv 

(uv-Yw)  +  2X4X7 


/u2~vv 

\uw“vw. 


—  + 


1  /u2-yv\  1X3  /u2-YvV 

2x9\uw-vw y  +  2x4x9 \uw-vw/ 


(46a) 


From  a  comparison  of  each  infinite  series  of  equation  (46a)  with  the 

series  expansion  for  (l-x2)  2  (formula  754  of  reference  7),  the 

following  equation  may  he  deduced: 
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1 


2X3X2X1  a3  ,  ^ 

"  TaT  ?  <‘V'T»> 

2*2*2 

* 


s®  ds 


(where  s  is  simply  an  integration  variable)  or 


jtcpn  =  2 


r(n-H)r  (|)an('TB-ug)n  C  \uv-vw  gZn 


H) 


^l-s2 


(46b) 


(The  integrals  of  equation  (46b)  may  be  expressed  as  incomplete 
beta  functions.) 


Along  the  boundary  U2  =  iv,  the  upper  limit  of  the  integra¬ 
tions  becomes  unity  and  the  values  of  the  integrals  are  given  by 
formula  482  of  reference  7  as 


r  HM§) 

2T(n+l) 
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In  this  special  case,  then,  cpa  reduces  to 


OO 

«<¥>a  =  "  ^  “  (vw-uw)n  =  rtf  (y) 


n=0 

as  is  required  along  the  curve  U2=uw. 

The  portion  of  the  x  component  of  the  perturbation  velocity 
associated  with  f(y)  may  be  obtained  by  partial  differentiation 
of  equation  (46b)  to  give 


^  _  _M_ 
bx  2Prt 


(-l)rT(n+l)r(j)an 


n=0 


H) 


Mn 


(u2-vw) 


1 

n~2 


(47) 


The  series  of  equation  (47)  is  a  function  only  of  vw.  Alterations 
of  the  function  f (y)  for  a  given  plan  form  cannot  change  the  factor 

When  an  arbitrary  strength  of  vorticity  is  allowed,  the  x  component 
of  the  perturbation  velocity  may  be  obtained  by  addition  of  equa¬ 
tions  (7a)  and  (47): 

r 


£>  _ 

he  "  2prf 


dv-du 


ab 


^(uv-u) (vw-v) 


Ua 

Prt 


hs) 

^Juw“u2 


^|vw-vi(u2)  - 


M 

2Ua 


n=0 


or 


Ua 
2P  it 


dv-du 


nr(n+l)r(|)an 

n-i 

(u2-Tw)  2 

rH)M” 

(48) 

(1-  ^)s(v„) 

(48a) 

^ab 


l\  (uy-u)  (vw-v)  ^uv_u2 
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TJcg 

where  g(vw)  represents  ^  times  the  braced  portion  of  equation  (48) 

The  function  g(v^)  depends  on  the  shape  of  the  wing-plan  boundaries 
(fig,  4(f)).  and  on  the  amount  of  vorticity  in  the  wake  of  the  sub- 

_ ,  .  ...  L  d^2\  ,  4  . 


sonic  trailing  edge.  The  quantity 
lines  of  constant  v™. 


g(vw)  is  constant  along 


If  solutions  that  satisfy  the  Kutta-Joukowski  condition  are 
desired,  the  constants  of  equation  (48)  may  he  evaluated.  Along  the 
subsonic  trailing  edge,  the  x  component  of  the  perturbation  velocity 
must  be  zero.  The  integral  of  equation  (48)  is  zero  along  the  boundary 
because  the  limits  of  integration  then  coincide.  The  Kutta-Joukowski 
condition  can  then  be  satisfied  only  if 


vw"vl(u2(vw. 


(-i)nr(n+i)r(|) 

rH) 


%  (u2-vw)n  2  (49) 


Equation  (49)  allows  all  the  coefficients  an  to  be  chosen.  Because 
the  first  member  is  finite  when  u2  =  vw,  ao  =  0.  The  evaluation  of 

of  the  other  coefficients  may  be  simplified  by  setting  p  =  j^u2-vw. 
Equation  (49)  then  becomes 

2Ua  ,| - 7 — 7 — i  \  (_1)nr(n+1)r(2)  an  ?  n 

—  — Tjppj —  ^  p2"-1  <««> 


fn  p2n-l 
Mn 


(49a) 


Equation  (49a)  is  a  power  series  in  p,  whose  nth  coefficient  an 
may  be  evaluated  in  terms  of  the  (2n-l)th  derivative  about  p  =  0, 
Differentiation  of  the  first  member  of  equation  (49a)  may  be 
accomplished  by  successive  applications  of  the  relation 


2  jjugjvy )-vw  d 

~du2(vv)  dvw 

dvw 


By  application  of  equation  (50)  and  L ’Hospital's  rule  to  equa¬ 
tion  (49a),  the  coefficients  aq  and  a2>  for  example,  have  been 
evaluated, 

dvi  du2 


=  -Ua 


du2  dvw 


lim  u2— »vw— »0  (51) 
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(• 


and 


The  other  coefficients  may  be  similarly  evaluated.  Once  the  coeffi¬ 
cients  ab,  &2)  ♦  •  •  ®n  have  been  determined,  all  the  perturbation- 
velocity  components  in  the  plane  of  the  wing  may  be  derived  by  means 
of  the  equations  previously  presented. 

A  special  case  of  equation  (5l)  gives  equation  (30) .  The  solu¬ 
tion  for  the  swept  trapezoidal  wing  may  thus  be  obtained  from  the 
general  equations  that  have  been  derived  without  the  a  priori  assump¬ 
tion  that  the  flow  is  conical. 

If  only  functions  derivable  from  the  x  component  of  the  pertur¬ 
bation  velocity  are  required  (such  as  pressure,  lift,  and  drag  coeffi¬ 
cients),  equation  (3),  (6),  (7),  or  (48)  may  be  directly  applied.  The 
solution  for  the > wing  of  figure  4(f)  that  satisfies  the  Kutta- 
Joukowski  condition  may  be  obtained  by  substituting  equation  (49)  in 
equation  (48)  to  yield 

Sce=u£L|  - —  (53) 

*  2f3*Jab  fuw-u  )(vw-v) 

For  the  special  case  of  a  wing  with  a  straight  leading  edge  (fig.  3), 
equations  (53)  and  (8)  give  equation  (12),  which  was  previously 
obtained  by  intuition. 


RECAPITULATION 

A  recapitulation  of  the  formulas  for  the  x  component  of  the 
perturbation-velocity  potential  of  a  single  thin  flat-plate  wing  tip 
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with  arbitrary  plan  boundaries  is  given  herein.  The  lift-distribution 
coefficient  Cp^L  of  thin  wing  tips  may  be  obtained  from  these 
formulas  by  the  relation 


(-2Cp) 


_  4  &£ 
"  U  ck 


(a)  In  regions  not  influenced  by  the  wing  tip. 


(from  equation  3(a)  and  fig.  1(a)). 

(b)  In  regions  influenced  by  subsonic  leading  or  trailing  edges 
when  no  vortex  sheet  exists  directly  behind  the  subsonic  trailing 
edge, 


(from  equation  (7a)  or  (48)  and  fig.  2(b)). 

(c)  In  regions  Influenced  by  subsonic  trailing  edgeB  when  an 
undetermined  amount  of  vorticity  is  allowed  in  the  vortex  sheet 


(from  equation  (48a)  and  fig.  4(f)). 
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(d)  In  regions  influenced  by  subsonic  trailing  edges  when  the 
Kutta-Joukowski  condition  is  satisfied  at  the  trailing  edge, 

r 


-/ab 


dv-du 

^|(uw-u)(vw-v) 


(from  equation  (53)  and  fig.  4(f)). 


Flight  Propulsion  Research  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Cleveland,  Ohio,  January  12,  1948. 
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APPENDIX  A 
SYMBOLS 

The  following  symbols  are  used  in  this  report: 

a0>  al>  a2>  •  •  •  ^  coefficients  of  power-series  expansion  of 

f(y)  =  a0  +  a^'  +  a2y2  .  .  . 


CP 

pressure  coefficient  on  top  wing  surface 

CP,L 

lift-distribution  coefficient  (  -2  Cp) 

f(y) 

perturbation-velocity  potential  on  top  surface 

of  diaphragm  in  wake  of  so-called  subsonic 
trailing  edge 

S  Ow) 

undetermined  function  of  vw 

kl 

constant  greater  than  0 

k2 

constant  whose  value  lies  between  0  and  1 

M 

free -stream  Mach  number 

m,  n 

integer  summation  indices 

P 

= 

^  u2"yw 

S 

plan-form  area 

u 

free-stream  velocity 

it 

3 

d-Pn)l 

oblique  coordinates  whose  axes  lie  parallel  to 

M 

v  "  2,6 

■(  t+Pn  )J 

Mach  lines 

x,  y. 

Z1 

Cartesian  coordinates 

t,  n. 

U 

a 

angle  of  attack 

M  .  \ 

P 

cotangent  of  free-stream  Mach  angle  (M M2-l j 

r 

gamma  function 
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X 


a 


<P 


slope  of  streamlines  near  z  =  0  plane  measured 
in  y  =  constant  planes  (diaphragm-surface 
slopes) 

slope  of  streamline  (on  wing  slopes)  with  respect 
to  z  =  0  plane  measured  in  y  =  constant 
planes 

perturbation-velocity  potential  on  top  wing 
surface  or  diaphragm 


Subscripts : 

I,  2,  .  .  • 

a 

b 

B 

D 

T 


w 

Examples : 

or 

xi 

CRD 

Sw(l+2) 

n 

r(n+l) 


% 


Si 


numbered  areas  or  curves 
portion  associated  with  f (y) 
portion  not  associated  with  f (y) 
bottom  (wing  or  diaphragm  surface) 
diaphragm 

top  (wing  or  diaphragm  surface) 
wing 

slope  on  top  wing  surface 

slope  of  diaphragm  in  plan-area  1 

perturbation- velocity  potential  on  top  surface 
of  diaphragm 

wing  area  1  plus  wing  area  2 
curve  v  =  v^(u) 
n  l 

portion  of  perturbation-velocity  potential  on 
wing  surface  associated  with  function  f(y) 

portion  of  diaphragm  slope  associated  with 
function  f(y) 

value  of  i  along  leading-edge  boundary  curve  1 
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APPENDIX  B 

DERIVATION  OF  EQUATION  (6) 

The  velocity  potential  at  point  (x,y)  for  the  wing  tip  of 
figure  2(a)  is  given  by  equation  (5)  (written  in  Cartesian  coordi¬ 
nates)  as 


The  area  a'b'd'c'  may  be  subdivided  into  the  areas  f  gkc ' ,  a'egf, 
and  eb'd'k.  The  area  fgkc'  is  the  same  as  the  area  abed  and 
the  relation  of  the  distances  in  the  denominator  is  the  same  for 
corresponding  points.  The  first  integral  of  equation  (B2)  may  there¬ 
fore  be  written 
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(B3) 

On  the  other  hand,  the  area  b'od'  may  be  expressed  as  the 
area  ghk  plus  the  area  eod'khg  minus  the  area  ebfdTk.  The  area 
ghk  is  the  same  as  the  area  bod  and  the  relation  of  the  distances 
in  the  denominator  is  the  same  for  corresponding  points.  The  second 
integral  of  equation  (B2)  may  therefore  be  written 
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The  first  two  integrals  of  equation  (B5)  are  surface  integrals  over  the  areas 
and  Sy  2>  respectively.  Integration  of  the  third  and  fourth  integrals  with  respect  to  \  over 
a  strip  of  width  dx  removes  one  integral  sign  and  replaces  d£  hy  dx.  The  width  of  the 
strip  in  the  fifth  integral  along  the  f  direction  may  he  determined  with  the  aid  of  the 
following  sketch  of  the  area  dd'k: 
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a  I  =  2(dx-dx2) 


dx2 


where 


is  the  slope  of  the  wing-boundary  curve  in  the 


vicinity  of  the  point  d.  Integration  of  the  fifth  integral  of 
equation  (B5)  with  respect  to  f  removes  one  integral  sign  and 
replaces  d|  with  the  fourth  member  of  equation  (B6) .  Equa¬ 
tion  (B5)  thus  becomes 
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If  equation  (B7)  ia  rewritten  in  oblique  coordinates,  equation  (6) 
results : 
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APPENDIX  C 

CALCULATION  OF  PRESSURE  COEFFICIENT  FOE  WING 

OF  FIGURE  3  BY  METHOD  OF  REFERENCE  1 

When  O-q  =  -Op  =  a,  equation  (5)  gives  for  the  velocity  potential 
on  the  top  surface  of  the  wing  shown  in  figure  3 


cp  -  ns 

^  Mn 


ruw 


U2(vw) 


/“'U 


2XJa 

m 


2Ua 

Mrt 


v 


du 

vw+kiu 


/~vw 


-kiu 


dv 

ftp 


du 


A/uw-u 

U2(vw)  » 

Al  [uw-«2(vw)]  [vw+kiu2(vw)] 


(kiUw+vw) 

+  -  tan 


klK'u2tYj] 

vw+kiu2(vw) 


Now 


Substitution  of  equation  (Cl)  into  equation  (C2)  yields 


£- 


(ci) 


(C2) 


Ua  , 

<lu2(vw) 

I 

A 

vw+kiu2(vw)  (kp+1)  . 

kl[uw-U2(vw)] 

1  "  dvV 

1 

uv-u2(vw)  h  tan  ^ 

vw+kiu2(vw) 

(C3) 
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The  pressure  coefficient  on  the  top  wing  surface  is  then 
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Figure  4*  -  Integration  boundaries  for  evaluating  velocity 
potential  in  regions  influenced  by  subsonic  trailing 
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(d) 

Figure  4*  -  Continued#  Integration  boundaries  for  evaluating 
velocity  potential  in  regions  influenced  by  subsonic  trailing 
edges# 
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Figure  4*  -  Concluded*  Integration  boundaries  for  evaluating 
velocity  potential  in  regions  influenced  by  subsonic  trailing 
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Figure  7,  -  Variation  of  y  component  of  perturbation  Telocity  in  plane  of  wing  near  tip 

of  swept  trapezoidal  wing. 


Pl^ 


Figure  8.  -  Variation  of  slopes  of  streamlines  in  plane  of  wing 
near  tip  of  swept  trapezoidal  wing. 


